ON THE SPEED OF APPROACH TO EQUILIBRIUM 
FOR A COLLISIONLESS GAS 



KAZUO AOKI AND FRANQOIS GOLSE 

Abstract. We investigate the speed of approach to Maxwellian equilibrium 
for a collisionless gas enclosed in a vessel whose wall are kept at a uniform, 
constant temperature, assuming diffuse reflection of gas molecules on the vessel 
wall. We establish lower bounds for potential decay rates assuming uniform L v 
bounds on the initial distribution function. We also obtain a decay estimate in 
the spherically symmetric case. We discuss with particular care the influence 
of low-speed particles on thermalization by the wall. 



In memory of Carlo Cercignani (1939-2010) 
1. Introduction 

Relaxation to equilibrium is a fundamental issue in statistical mechanics. In the 
kinetic theory of gases, convergence to equilibrium is formulated in terms of the 
asymptotic behavior of solutions of the Boltzmann equation in the long time limit. 
In Cercignani formulated a quite remarkable conjecture, asserting the domina- 
tion of Boltzmann's relative H functional by a constant multiple of Boltzmann's 
entropy production functional. If true, Cercignani's conjecture would imply expo- 
nential relaxation to equilibrium for space homogeneous solutions of the Boltzmann 
equation, and would have important consequences even for space inhomogeneous 
solutions. This conjecture was subsequently disproved in a series of papers by 
Bobylev, Cercignani and Wennberg [2j |3j [14]. And yet, in the words of Villani, 
"Cercignani's conjecture is sometimes true, and almost always true" [12] . Slightly 
later, Desvillettes and Villani [7] studied the speed of convergence to equilibrium 
for space inhomogeneous solutions of the Boltzmann equation, under rather strin- 
gent conditions of regularity, decay at large velocities and positivity, but without 
assuming that the initial data is close to cquilbrium. They found that the distance 
between the distribution function at time t and the final equilibrium state decays 
like 0(t~ m ) for each m > as t — > +oo. Their result makes critical use of the Boltz- 
mann collision integral and entropy production, and of their combination with the 
free transport operator (the left-hand side of the Boltzmann equation.) In an ear- 
lier paper [5], Desvillettes had proved that Boltzmann's entropy production can be 
regarded as a distance between the distribution function and the set of Maxwellian 
states. On the other hand, Guo proved in [5] that, for initial data close enough to a 
uniform Maxwellian, solutions of the Boltzmann equation converge to equilibrium 
at exponential speed. 
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In [7] , Desvillettes and Villani treated the cases where the spatial domain is the 
flat torus (equivalently, the Euclidian space with periodic initial data) , or a smooth 
bounded domain with either the specular or bounce-back reflection conditions at the 
boundary. Therefore, convergence to equilibrium is driven by the collision process 
alone, and the role of the boundary condition is rather indirect. In the long time 
limit, the solution of the Boltzmann equation converges to a uniform Maxwellian 
state whose density and temperature are given in terms of the initial data by the 
conservation of total mass and total energy. The boundary condition is involved 
only through its compatibility with the conservation laws of total mass and energy. 
In a more recent work [13] , Villani extended their result to the case where the gas is 
enclosed in a vessel whose walls are maintained at a constant temperature, modeled 
with an accommodation or diffuse reflection condition. In that case, mass is the 
only quantity that is conserved under the dynamics of the Boltzmann equation, so 
that the final equilibrium state is the centered, uniform Maxwellian with the wall 
temperature and same total mass as the initial data. In that case, the boundary 
condition has a direct influence on the relaxation to equilibrium. Nevertheless, the 
combined effect of entropy production by the collision process and of free streaming, 
an example of what is termed "hypocoercivity" , is the driving mechanism used by 
Villani in his proof of relaxation to equilibrium [13] . 

However, it seems obvious from physical consideration that the accommodation 
or diffuse reflection condition on a wall maintained at some constant temperature 
is enough to drive the gas to an equilibrium state, and that the role of the collision 
process is not as important as in the periodic, or specular or bounce-back reflection 
cases — which are in any case less realistic. 

In order to establish the role of Boltzmann's collision integral in the rate of 
relaxation to equilibrium, we consider the case of a collisionless gas enclosed in a 
vessel whose surface is maintained at a constant temperature, and investigate the 
speed of approach to equilibrium for such a system. The present paper provides 
a theoretical confirmation of the numerical results obtained by the authors in [10] 
in collaboration with Tsuji. Basically, while the collision process is not necessary 
in the convergence to equilibrium, it greatly influences the speed of relaxation to 
equilibrium. In the same spirit, Desvillettes and Salvarani [6] have investigated 
the speed of relaxation to equilibrium in the case of linear collisional models where 
the collision frequency is not uniformly bounded away from 0, by adapting the 
hypocoercivity method. 

We dedicate this piece of work to the memory of our teacher, friend and colleague 
Carlo Cercignani, in recognition of his leading role in the mathematical theory of 
the Boltzmann equation. 



Consider a collisionless gas enclosed in a container materialized by a bounded 
domain ft in R w whose boundary d£l is maintained at a constant temperature 
6 W > 0. Its distribution function F = F(t, x, v) satisfies 



2. Statement of the problem 



( d t F + v ■ W X F 



0, 



(x,v) e n x r n , 




F(t,x,v')v'-n x K(x,v,dv') , (i,»)6L, (1) 



= F' 
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The notation n x designates the unit outward normal at the point x £ dCl, and 
we have used the notation 

T+ := {(x, v) £ dfl x K N | v ■ n x > 0} , 

T_ := {(x, v) £ dVL x K N | v ■ n x < 0} . 

The boundary condition is given by the measure-valued kernel K(x,v,dv') > 
satisfying the assumptions 

(i) for each x £ 9f2, one has 

K(x, v, dv')dv = dv 

v-n x <0 

(ii) for each x £ dSl, one has 

M(i t a t e)(v)\v-n x \ = / M(ifl t e){v')v'-n x K(x,v,dv') , v-n x <0, 

Jv'-n„>0 

iff 9 = 9 W . Throughout this paper, we use the following notation for the Maxwellian 
distribution: 

M M) {v) = ^J^e-^ 

(Henceforth, we denote for simplicity M w = M.n t o t e w )-) With these notations, the 
case of diffuse reflection corresponds with 



K(x,v,dv'):= ^o,e w) (v)\v-n x \dv' 



M(i t a t e w )(u)\u ■ n x \du 

u-n x <0 



The case of an accommodation boundary condition (i.e. the so-called Maxwell-type 
condition) corresponds with 

K (x, v, dv') := (1 — a(x))5(v — v + 2(v ■ n x )n x ) 
M( lfi M m )(v)\v ■ n x \dv' 



a(x) 



■M.(i^ 8w )(u)\u ■ n x \du 

u-n x <0 



In the case of diffuse reflection, Arkeryd and Nouri [T] have proved that the 
solution of (UJ satisfies 

F(t+ •,-,•) -t f f uinj^j^, ;„ rl/Tn T-l ^ O ~ n N \ 



mi 



ff F m dxdv in L x ([0,T] x Q x 
JJnxR N 



as t — > +oo. 

In the present paper, we study in detail the following 

Problem: Is there a decay rate E(t) — > as t —> +oo, for which an asymptotic 
estimate of the form 



SlxR" 



= 0(E(t)) (2) 



holds for each solution F of (p} uniformly as F m runs through certain classes of 
bounded subsets of L 1 (R N x R w )? 
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3. Negative results and lower bounds on E 



We first consider the problem of finding a decay rate that is uniform as F m runs 
through all bounded subsets of L 1 (R 7V x R ff ). 



Proposition 3.1. There does not exist any decay rate function E : R^ 
such that E(t) — ¥ as t — > +oo and 



F(t + -,v) 



F m dxdv 



n x r 



LUlO.TlxfixR") 



R4 



(3) 



< E{t)\\F v 



for each t > 0. 



Since the initial boundary value problem (p} is linear, if there is a long time 
decay as in ^ uniformly as the initial data F m runs through bounded subsets of 
L , an inequality of the form Q must hold. 

Proof. Assume that the origin € il — if not, one can apply a translation in the 
x- variable without changing the nature of the problem — and set R — \ dist(0, dil). 
We henceforth denote by B the closed unit ball of R N , and by 1a the indicator 
function of the set A. 
For < e < 1, set 

F m (x,v) = ^1 £ b(.t)1 e b(«) 
(with eB :— {ez \ z G B}) so that 

//^<*>»>«*««=|B|>. 

Call t(x,v) := mf{t > | x — tv ^ O}; obviously the solution F of ([T|) satisfies 
F > $, where $ is the solution of the same free transport equation, with the same 
initial data, but with absorbing boundary condition: 



f d t <f> + v-V x $ = 0, (x,v) e n x R^, 
$(i,a;,w) = 0, e T_ , 



(4) 



t=o 



Solving for $ along characteristics, one finds that 



Now 



F(t 



^(l,0.9 m ) 



F m dxdv 



flxTL N 



> 



F(t + 

f [ [ (F(t + s,x,v)- 
Jn Jr n \ 



IBI 2 

)-yr^(w 



n| 

]BJ! 
mi 



L 1 ([0,T]xQxR N ) 



L 1 ([0,T]x[JxR w ) 



\ + 

A / I(i.o,e„)('f) I dvdxds . 
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Henceforth, if A C R 3 and u <E R 3 , we denote A + u := {z + u \ z € A] . Since the 
function zh>z+ = max(z, 0) is nondecreasing, one has 



SlxR" 



> 



Jo In In* (*( < + s ' a; ' u ) ~ ^ M (^ w )(v)j dvdxds 

mi "W leB+(t+s)« (^)leB (") lo<t+s<r(x,v) 

|B| 2 \ + 

— \n[ M{1 - ' e,jo) ( v y dvdxds 



rT , , / j 



- ~2]v 1 £B + (t+ ;s ) 1 ;(^)leB(v)lo<(t+ S )| 1 ;|<i?-e 

Jo Jr" V £ 



IBI 



"•^ / '(i,o,e„)(' !; ) ) dvdxds 



\n\ 

since r(x, > 2i? — e for each x e eB. Now 

/ 1 |B| 2 

( ~^2N 1 eB+(t+s)t)(a;)l£B(w)lo<(t+s)|i;|<i?-e 

1 / e 27V |B | 2 V 

= ^■'i-e-B + (t+s)v(x)le-B(v)lo<(t+s)\v\<R-e ( 1 i^i ^(l,0,e w ) («) ) 

1 . . . . / £ 2N |B| 2 \ + 

l f l -^N 1 EB+(t+s)v(x)leB{v)lo<(t+ s )\v\<R~edvdxds 

Jo Jq Jr n e 

' B l / / -M 1 e-B{v)l <(t+s)\v\<R-edvds 
JO JR N e 



Then 



— 1^1 / / _ iv : '-0<|u|<mm(e,(_R-e)/t+s)^rfs 
JO JR N £ 



|B| rimin^l^U, 
>|B|4min/> (i? ~ £)JV 



Therefore, if there is a uniform decay rate £ as in ([3]), it must satisfy 

Here £ only has to satisfy < £ < R. If one lets £ — > + , this contradicts the 
existence of a profile E in ([3]) such that E(t) — » as i — > +oo. □ 

The negative result in Proposition l3.1l could have been easily anticipated: as e — » 
0, our choice of F m converges to |B| 2 <5 K= o^=o in the weak topology of measures. 
Then particles distributed under S x= qS v= o do not move, hence cannot reach the 
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boundary and therefore cannot be thermalized. Thus, in some sense, the problem 
of finding a decay rate as in ([3j> is trivial if stated in this way. 

Before going further, let us comment on the kind of inequality considered in 
Proposition 13.11 In order to measure the speed of approach to equilibrium for the 
problem ([1]) , it might be more natural to consider a pointwise in time inequality of 
the type 



S!xR f 



<E(t)\\F m \\ Li{nNxRN) (5) 

instead of ([3]). Observe that the l.h.s. of ([3]) is the time-average of the l.h.s. of 
([5]), so that (0 is a stronger statement than ([3]) — as a matter of fact, ([5]) entails 
([3]) with E(t) changed into TE(t). Therefore, negative results obtained on time- 
averaged estimates like ([3]) are stronger than their analogues for the corresponding 
pointwise estimate similar to (JSJ. This is why we sought to confirm numerically 
pointwise in time bounds like ([5]) in our previous paper [10] . and not time- averaged 
estimates as in the work of Arkeryd-Nouri pQ. Throughout the present section, 
we discuss negative results on the speed of approach to equilibrium and therefore 
consider time-averaged decay estimates. 

A slightly less trivial issue is to investigate the existence of such a decay rate 
under the additional assumption ruling out the possibility of a complete concen- 
tration as above. This could be done by controling not only the L 1 norm of the 
initial data, but (for instance) its entropy, or relative entropy with respect to the 
final state. For instance, one could look instead at a bound of the form 



F(t+;;.)- M ^' II F-.inlr 



SlxR" 



i 1 ([0,T]xnxR' v ) 



<E(t) [[ 



(G) 



F in {x,v)\\nF m {x,v)\dxdv . 



Proposition 3.2. If there is a time decay rate E such that the solution of (QJ) 
satisfies (0|), then it must satisfy 

C 

E(t) > as t — >• +oo . 

mt 

Proof. For the same initial data as before 

F m (x 7 v)\\nF m {x,v)\dxdv = J J ^\ eB {x)l eB {v) In dxdv 

= 2iV|B| 2 |lne|. 

Proceeding as before to obtain a lower bound for the l.h.s. of (|6]), one arrives at 
the inequality 

TIB! (l- , f N ^ M/ , ) (mm(l, (?~ e \ X) < 2A^|B| 2 |lne|£;m. 
Choosing £o > small enough, this inequality takes the form 

c(mm(l,^\\ <\\ne\E{t) 
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for all t > T and all e G (0,£o)- Picking then e = 1/t with i > 1 so that in 
particular t > max(T, 1/eo) leads to 



C 

E(t) > — as i ^ +00 . 



□ 



The constraint on the decay rate E obtained in Proposition l3.2l rcmams the same 
if one replaces the r.h.s. in the bound ((6]) with either 



E{t) // F m (x,v)(l + \v\ 2 + \lnF m (x,v)\)dxdv 
E(t) [[ ( F in (x,v)ln(—[— ) -F m +pM {lW ) (x,v)dxdv. 



where 



p=-ll F-dxdv. 



Indeed, compared to the r.h.s. in the bound ©, both these quantities involve in 
addition the total energy integral 

\v\ 2 F m (x,v)dxdv . 

fixR N 

Since the inequality © is tested with velocity distribution functions concentrated 
near v — 0, this additional term is negligible and its presence docs not change the 
final result. 

Of course, prescribing an entropy bound is not the only possibility for avoiding 
concentrations in the initial data F m . For instance, one could ask for a control of 
F m in some L p -norm with 1 < p < 00. In other words, one could investigate the 
possibility of an inequality of the form 



F(t + .,;.)- M ^> jj F-.lnlr 



SlxR" 



i 1 ([0,T]xfixR w ) (7) 



< E (t)\\ Ft 



Lp(OxR n ) 



Proposition 3.3. // there is a time decay rate E such that the solution of ([JJ) 
satisfies then it must satisfy 

C 

E {t) > AT ■ „ nl M aS t +°° • 

v ' — fN min(l ,2 / p' ) 

Proof. With the same initial data as before 

l B | 2/p 

\\F m \\Lp(nxR») = £2N/p , , 

with p' = < 00 the dual exponent of p, so that, estimating the l.h.s. by the 
same argument as before, one arrives at the inequality 



2AT| 



B| 2 \ ( . ( (R-e) \Y |B| 2/p 



s 
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Assume that p > 2. Since the growth of the l.h.s. in 1/e is smaller than that of the 
r.h.s., there is no point in letting e — > 0. One therefore chooses e > small enough 
but fixed, so that e < R/2 and e 2N \B\ 2 < \\VL\(2irO w ) N / 2 , so that 

In that case, one sees that there exists C = C(T, R, e) > such that 

C 

E{t) >-pj as t > 1 . 

If 1 < p < 2, by picking e > small enough, the inequality above takes the 
form 

1\\ N \B\ 2 /p 



for all i > and all e € (0, £o). Picking e = l/t with t 3> 1 so that in particular 
t > max(T, 1/so), one concludes from this last estimate that 

C 

□ 



At this point, it is worth comparing the result in Proposition 13.31 for the col- 
lisionless gas driven to equilbrium by the boundary, and the following important 
result, obtained by Ukai, Point and Ghidouche for the Boltzmann equation 
linearized at a uniform Maxwellian equilibrium. 

Consider the Boltzmann equation linearized at the uniform Maxwellian state 
M = A4 (1,0,1) j assuming hard sphere collisions, in the case where the spatial domain 
is the periodic box with unit side T 3 . In other words, consider the Cauchy problem 

f (d t +vV x )g+£ M g = 0, (t, x, v) e R* + x T 3 x R 3 , 

The linearized collision operator is given by 

C M <t>(v) := // O(v') +</>«) ~ ( t ) { v ) ~ ~ «*) • u\M(v*)dv*du 

J JR 3 xS 2 

with 

{V 1 = V — (V — • LJLU 

V* = + (v — V*) ■ LdUl 

where cu runs through the unit sphere S 2 . It is understood that the notation 
^M<^(i) x, v) designates (CM<fi(t, x, -)){v). Define 

n 5 m = (g m ) + (vg m ) ■ v + ((1M 2 - l)g in )\{\v\ 2 - 3) 

with the notation 

((f)) = / <j){v)M{v)dv . 

JR 3 
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Theorem 3.4. There exist 7, > 0, and, for each"/ €]0, 7*[ } a positive constant 
C-y such that, for each g m £ L 2 (T 3 x R 3 ; Mdvdx), the solution g of (Bty satisfies 
the inequality 



<?(*,-,•) 



T 3 



Ug m dx 



< C 7 e 7t ||5 m ||L2( T 3 xR 3. M ^ da ,) 



L 2 (T 3 xTL 3 :Mdvdx) 

for each t > 0. 

This result expresses that there is a spectral gap for the linearized Boltzmann 
(unbounded) operator 

Lg(x, v) = v ■ W x g(x, v) + C M g(x, v) 
on L 2 (T 3 x R 3 ; Mdvdx) with domain 

D(L) = {<j)e L 2 (T 3 xR 3 ;Mdvdx)\v-y x (t) and (l + \v\)<p e L 2 (T 3 x R 3 ; Mdvdx)} . 

In other words, the spectrum of L is included in a domain of the form 

spec(L) C {0} U {z e C I Re(z) > 9} 

for some 6 > (the spectral gap). The eigenvalue corresponds to the five- 
dimensional nullspace of L, given by 

kerL = {a + b ■ v + c\{\v\ 2 - 3) | a,c € R and b e R 3 } . 

One reason explaining the spectral gap as above is the fact that Cm is an (un- 
bounded) Fredholm operator in the v variable, of the form 

£ M <P(v) = v{\v\)<t>(v) - K<t> 

where ^(|f|) is a multiplication operator by the collision frequency satisfying the 
inequality u*{l + \v\) < ^(|w|) < v*{l + \ v\) for some positive constants and v* . 

Proposition 13.31 for p = 2 rules out the possibility of a similar spectral gap 
property for the problem ([1]). 

Indeed, assume for simplicity that 6 W = 1 in the problem ([TJ. Then by the 
Cauchy-Schwarz inequality 

F(t,;-) 1 



> 



M 



— // F m dxdv 



L 2 (fixR3 -Mdvdx) 



M 

W\ 



F m dxdv 



OxR 3 



so that a spectral gap inequality of the form 

F(t,;-) 1 



M 



— / / F m dxdv 



OxR 3 



M 



would entail 



F(t 



•)-- 



F m dxdv 



fixR 3 



F 



L 2 (nxR 3 -Mdvdx) 
L 2 {nxR 3 ;Mdvdx) 

L 1 ( [0 , T] x a x R 3 ; dv dx ds ) 
L 2 (nxR 3 ;Mdvdx) 



(9) 



M 



10 



K. AOKI AND F. GOLSE 



Now, for 

F in (x,v) = ^l £B (x)l eB (v) 



one has, assuming that < e < 1 
pi 



M 



< ((2n) N e) 1/4 \\F m \\ L 2(n x - R 3. dvdx \ 

L 2 {SlxR 3 ;Mdvdx) 



In other words, the existence of a spectral gap for the problem (TjQ) in il x R 3 
i.e. the inequality @ would entail the inequality 



S!xR 3 



L 1 ([0,T] X SI x R 3 -dvdxds) 



which stands in contradiction with Proposition 13. 31 

4. Algebraic decay 

After the negative results obtained in the previous section, we now consider the 
problem of establishing decay estimates compatible with the results of section [31 

We shall do so on a special class of boundary value problems, satisfying the 
following assumptions: 

(a) the spatial domain 51 is a ball of the Euclidian space R 3 ; 

(b) the distribution function of the collisionless gas at t = is a radial function of 
both the space and velocity variables; 

(c) the boundary condition is a diffuse reflection (complete accommodation), with 
a constant boundary temperature (no temperature gradient at the boundary of the 
vessel.) 

None of these assumptions excludes the mechanism producing the negative re- 
sults in section |3l Indeed, all these results are based on the choice of an initial 
distribution function that is concentrated, in the velocity variable, near v = 0, and 
in the space variable, near a point that is away from the boundary. This is obviously 
possible with the symmetries above, by choosing the initial data concentrated, in 
the space variable, near the center of the ball, and in the velocity variable, near 
v = as in section [3] Particles so distributed at time t = will not reach the 
boundary of the spatial domain until after very late, and, meanwhile, the exact 
nature of the boundary condition is immaterial for such particles. 

Therefore, we set the spatial domain to be ft = B(0, 1), the unit ball of R 3 . 
The temperature on dQ is set to be 8 W = 1. Consider the initial boundary value 
problem 

(d t + v ■ V x )f = 0, (t,x,v) G R+ x n x R 3 , 

/(t ) i,«) = ^M(ti) ) vn x <0, x&dn, (10) 
\ M ) + 

with the notation 

(</>)+= / <j>(v)v ■ n x dv , 

and we recall the notation M — A^i^i) for the Maxwellian distribution with 
density 1, zero bulk velocity and temperature 1. 
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The main result in this section is the following. 

Theorem 4.1. Assume that < f" l (x,v) < CM(v) a.e. in (x,v) e fixR 3 . Then 
the solution of U 0\) satisfies 



f(t,x,v) x 



so that 

/(*,'»•) i 



f m {\x\,\v\)dxdv 

!)xR 3 



1 \ „ f 1 



M l°l 
for all 1 < p < oo . 



/* B (|a:|,|t)|)da:dt; 

!!xR 3 



= o 

LP(axR 3 ;Md?;d2;) 



1 



£min(l,3/p) 



The reader might want to compare the upper bound in this theorem with the 
lower bound in Proposition 13.31 this is done by taking p — oo in the estimate ([7]) 
and in Proposition 13.31 and p = 1 in the theorem above. The difference in the 
exponents observed suggests that one of these results at least is not sharp. The 
numerical investigations in |10) suggests that it is the lower estimate in Proposition 
13.31 that is sharp, while the upper bound in the theorem above is not. 

Our main task in the proof of this theorem is to use the symmetries in the problem 
(|10[) to reduce it to a renewal integral equation for a scalar unknown quantity, and 
to use classical results for such equations. 

The role of symmetries in this problem is summarized in the following statement. 

Lemma 4.2. If f is the solution of the initial boundary value problem U0\) . its 
outgoing flux at the boundary 



fi := f(t, x, v)v ■ n x dv 

J v-n x >0 

is independent of the position variable x as x runs through dfl, and a function of 
t > only. 

Proof. If i?€ C>3(R) is any orthogonal matrix, the function x, v) := f(t, Rx, Rv) 
obviously satisfies /fl| t=0 = /| t=0 , and (f B ) + (t,x) = (f) + (t,Rx), while 

(d t + v V x )f R (t, x, v) = d t f(t, Rx, Rv)+v {R T V x f){t, Rx, Rv) 
= d t f{t, Rx, Rv) + Rv ■ V x f(t, Rx, Rv) = . 

By uniqueness of the solution of ([Toll , we conclude that fu = f for each R E 0$ (R) . 
Thus 

(f)+(t,x)={f R ) + (t,x) = {f)+(t,Rx) 
for each (t,x) £ R + x d£l and R E (^(R): hence /j, := (/)+ is independent of x as 
x runs through dQ. □ 

Proof of Theorem \4-l\ With Lemma 14.21 we easily arrive at an integral equation 
satisfied by fi. Pick x £ dil, and set 8 = (n x ,v); the backward exit time 

t x v := inf{t > I x — tv ^ fi} 

is given by the explicit formula 

2 cos 9 



T 



x,v 
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as can be seen on Figure 1. 

By the method of characteristics, for each x € dVL and v € R 3 such that v-n x > 0, 
one has 

f(t, X, V) = l t>Tmv ^-fi(t - T X , V ) + h<r x ,J in (\x ~ tv\, \V\) . 
\ M ) + 

Hence 

H(i) = / f(t,x,v)v ■ n x dv 

Jv-n x >0 

= TT7T / n{t - T x , v )l t>Tx v M(v)v ■ n x dv 
J v . nx>0 

+ / f in (\x-tv\,\v\)l t<rxv vn x dv 

J v-n x >0 

Introducing spherical coordinates in v, with n x being the polar direction and 9 the 
colatitude, we see on Figure 1 that, for each x € dfl, one has = 2cos0, so 

that 

/ fj,(t - T x v )l t>Tx v M(v)v ■ n x dv 

Jv-n x >0 

r°° r*/ 2 / 2cos#\ 
= 2ttJ J n^t — J l t{v>2cos e M{\v\)\v\ cos 6\v\ 2 sin 6d8d\v\, 



abusing the notation M(|u|) to designate the radial function v h-> M(v). 
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Substituting y — cos 9, this integral is expressed as 

cOO pi 



JO 



fi (t - ^ l tH>2y M(\v\)\v\ 3 ydyd\v\ 



2/M 

M(t-r)lt> T M(|w|)|t;| 5 T(iTdl«| 



^0 

oo pt 



H(t-T)l T<2/lvl M(\v\)\v\ 5 TdTd\ V \ 

JO 

= \f\{t-r) (j*'* M (M)M 5 dM^) rdr . 

We henceforth define 

K(r) ■■-^J o 2/T M(\v\)\vfd\v\. 

The function K is continuous on Rl, satisfies K(t) > for each r > 0, and 

g 

X(t) as r ->• + , if (r) ~ ^— ^ as r ->• +oo . 
Therefore ju satisfies 

H(t) = / K(T)n(t - r)dr + S(t) 
Jo 

where the source term is 



S(t) := / f m (\x - tv\, \v\)l t<T ^v • n x dv . 

Jv-n w >0 

This is a renewal equation, to which we shall apply the results of [8] . 
Next we consider the source term. Define 

$(M):=supe S s, en / in (|:i:|>|); 

assuming that < f m £ L°°(Cl x R 3 ) implies that $ £ L°°(R + ). Then, 



S(t) < / $(|w|)l t<T . , v v-n x dv 

J v-n x >0 

P Tr/2 



/•OO / PIT / A \ 

= 2tt y $(|w|)|w| 3 ly li|„|<2c OS e cos 6 sin 6W0 j d| 
= 2tt^ H\v\)\v\ 3 ^ l tlvl<2y yd y ^j d\v\ 



w $(\v\)\v\*{l-tf\v\ 2 )+d\v\ 



< 7r||$||jroo / \v\ d d\v 



2/t 



Apply Theorem 4 of [8]: with n = 3, one has 



/>oo 

/ t m K(t)dt < +oo for m = 0, 1, 2, 3 
Jo 
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while 

as t — > +00. Then 



t I . 4 \t 



fj,(t) -> fl 



S{r)dT 



OO poo 



poo 

/ T K(r)dT 
Jo 



as t — > +00, and 

- Moo = O ^ 

Now let 5 := f/M, where / is the solution of (|T0l) ; g satisfies 

(d t + w • V x )5 = 0, (t,x,v) e R+ x n x R 3 , 

<7(i, x, u) = v / 27r/i(i) , w • n-a < , x € 5f2 , (1 

=0 = / in (M>|)/M(V). 
The function g is given by 

f ,;,l flxl lull 

g(t, x, v) = V2^v(t - r x>v )l t>TxiV + ^gj ' ; l t <r x ,„ , 

so that 

X, U) - V^Moo = \/27T (/x(i - r X)V ) - /loo) lt>r x ,„ 

7 in (M,M) 



M(u) 



V^TT^oo ) lt<T x ,„ = I + II . 



Now Ta;,^ < 2/|u|, so that 



l t <r.,„M(t;)ii« < f M(v)dv ~ 
Hence 

1 



On the other hand 



\II\ p M{v)dv = 01-). (1 



\I\PM(v)dv < / ^_M(i,)d V 



-I Tti M Up M(v)dv+ f — 2 ° M(v)dv 

l\v\<e {{t ~ —\)+ + l) P J\v\>e [{t ~ + l) p 



c_ 



< C I M(v)dv + — — ^ — r - rT7 / M(v)dv 



(-11 

< C'e 3 + -j C 



for some other constant C > 0, provided that t > 4/e. Choosing z = 4/t < 1, 
arrive at 

43 (~ii ni 

|J|"Af (»)d« < —5- + — . 

/ R JV t t F 
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Putting together (|T3 |) - (fT4| entails the first estimate in the theorem. The second 
follows from the first and interpolation with the maximum principle for the solution 
g of CC2). 

Finally, the fact that 



27^=1// f m (\x\,\v\)dxdv 
follows from the conservation of mass for the solution of (fTUl). □ 



Notice that the decay requirements on the renewal kernel in Feller's Theorem 4 [8] 
forbid extending our proof to space dimensions less than 3. In the one-dimensional 
case, Yu obtained in [TS] an algebraic rate of convergence for the outgoing flux 
(specifically, of order t -1 / 10 ), by probabilistic arguments significantly more involved 
than the above proof. 

5. The monokinetic case: exponential decay 

The discussion in section [3] shows the existence of too many slow particles is 
the reason for the slow return to equilibrium in the case of a collisionless gas in a 
container with constant wall temperature. In the present section, we confirm this 
effect by proving that monokinetic populations of particles in the same conditions 
converge to equilibrium exponentially fast. 

Radiative transfer provides a natural example of a monokinetic transport process. 
Consider the following (academic) problem: at time 0, a source emits radiation in 
an empty container surrounded by infinitely opaque material. Photons propagate 
at the speed of light, and are instantaneously thermalized on the wall. Denoting 
I = I(t, x, uj, v) the radiative intensity at time t, position x, in the direction uj and 
with frequency v, one arrives at the following initial-boundary value problem: 

\d t i + uj ■ v x i = o , x e o , M = l , v > o , 

I(t,x,w,v) = B„(6i(t,x)) , ieffl,W'i! I <0,i'>0, (15) 

l\ = I in . 
\t=o 

Here c > is the speed of light in vacuum, f2 is a bounded open set in R 3 with C 1 
boundary, n x is the outward unit normal at a; € dfl. The notation B u {0) designates 
Planck's function for the radiative intensity emitted by a black body at temperature 
9, i.e. 

r (0) - 2hv " 1 

where h and k are respectively the Planck and the Boltzmann constants. The 
temperature 9i is defined by 

I(t, x, uj, v){ui ■ n x ) + du)di> 

S-XR+ 

B y {9i{t,x)){uj ■ n x ) + dujdv =: a9i{t,x) A 

/S 2 xR+ 

where a is the Stefan-Boltzmann constant 

27rfc 4 f°° x 3 dx 2n 5 k 4 



2 h 3 ./o e x - 1 15c 2 h 3 
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By definition of 9i, the net flux of energy at each point of the wall is 
/ / I(t, x, LO, v)lo ■ n x duodv = , 

J J S 2 xR + 

so that the total energy in the container is conserved: for each t > 0, one has 
I(t,x,u>,v)dxdujdv = / // I m (x,Lo,v)dxdLodv . 

fixS 2 xR + J J JnxS 2 xR + 

Therefore, one expects that, as t — > +oo, the radiation inside the container ther- 
malizes so that 

I(t,x,uj,v) — > B^^oo) as t — y -f-oo , 

with 



/■OO 

Jo 



SlxS 2 xR 



I m {x, LO, v)dxdiodv . 



+ 



In this section, we seek the rate of convergence to equilibrium for the problem 
above. 

Theorem 5.1. Assume that I m is smooth and satisfies supp(/™) C 51 x R + to- 
gether with the bound < I m (x, v) < £>„(0) with Q > 0, for all x € and v > 0. 
Then, the solution I of the initial boundary value problem \15)) satisfies 

I(t,x,uj,p)-B u {9 0O ) = 0(e- Q '*) 
as t -> +oo, uniformly on 57 x S 1 x R + , for each a' < a, where Qoq > is given by 

I m (\x\ 2 ,v)dxdv 



00 a\fl 
while a is 



SlxS 2 



a = inf |-3fi(z) | z G C* such that | J e sz sds = l| . 

Proof. The problem is simplified by considering the radiative intensity integrated 
over frequencies, the grey intensity u defined as 



u(t,x,io) := / l{t,x,Lo,v)dv . 
Jo 

The formula defining (9/ reads 

u(t, x, lo)(lo ■ n x )+dio = a9i{t,x) A 



s- 



while 



/ B y {6 I )dv = -9j. 
Jo n 

Therefore, choosing for simplicity units of time and space so that c — 1, the grey 
intensity u satisfies 

d t u + lo ■ V x u = , x e , \lo\ — 1 , 

u(t,x,Lo) = — / u(t,x,Lo')(Lo' ■ n x ) + duj' x <E dCl , lo ■ n x < , , „, 

n Js 2 U"J 

oo 

I in dv . 
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Henceforth, we assume the same geometric setting as in sectionS) = B(0, 1) C 
R 3 , and I m is radial in x and isotropic, meaning that I m = 7*™(|ar| 2 , v). In that 
case, reasoning as in Lemma 14.21 shows that 

/(*) := i / u(t,x,ui')(u;' -n x ) + du;' 
Js 2 

is independent of the position x € dCl. 

Proceeding as in section |4l we derive a renewal integral equation for /. Indeed, 
integrating the transport equation in (|16[) by the method of characteristics, we see 
that 

u(t,x,u) = U<r a , w u in (\x - tuj\ 2 ) + l t>r ^/(i - r XiU ) (17) 
whenever x £ d£l with oj ■ n x > 0, so that 

/(*) = b i u(t,x,w)(u-n x ) + du) 
Js 2 

lt<r x u , um (\ x ~ tui\ 2 )(u} ■ n x )+du) 

s 2 

+ b 1-t>T KU f(t — T x ,u){u •n x )+dio 
Js 2 

= S(t)+ [ K(s)f{t-s)ds. 
Jo 

Here 



so that 



while 



K(s)f(t-s)ds = ± / l t>Tx u f(t - t x>u ){u) ■ n x ) + duj 
Js 2 

«r/2 

= b lt>2 cos e f{t — 2 cos 9) cos 8 ■ 2tt sin 9 d9 
Jo 

= / U> s f(t - s)±sds 
Jo 



t 

K(s)f(t-s)ds 



K(s) = \s\ s<2 , 

S(t):=i / l t <T x y n {\x-tu\ 2 ){u-n x ) + du 
Js 2 

= b I 1 *<^ ^ m (! + t2 - 2tx ■ ' n *)+ duj 

Js 2 

,w/2 

b / lt<2cos0U m (l + t 2 -2t cos 9) cos 9 -2tt sin 9d9 
n Jo 

= \\ l t < s u m (l+t 2 -ts)sds 
Jo 



2 



U in (l-tT)K(t + T)dT. 



(I 
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Since u m £ L°°(Q), and K is bounded with support in [0, 2], so is S and one has 

/ t' m K(t)dt < +oo and S(t) = 0{r m ) as t -> +oo for each m > 1 . 
Jo 

Applying Theorem 4 of [8] already implies that 



S(t)ttt 

/(*) -> /oo := ^755 as f +00 . 

£i\:(i)<ft 



and that /(t) — is rapidly decaying as t — > +00: 

!/(*) - fool = 0{t~ m ) as t -> +00 for each m > 1 . 

We improve this estimate and obtain exponential convergence, as follows. Given 
4> € L 2 (R+), we denote its Laplace transform by 

/■OO 

4>(z) := / e- zt 4>{t)dt. 
Jo 

We recall that the Laplace transform </> 1— > is an isomorphism from L 2 (R + ) to the 
Hardy space H 2 (Ho), where IT Q :— {z G C 1 > a}. We recall that an element 
of W 2 (H a ) is a holomorphic function tp = ip(z) on the half-plane II a that satisfies 

sup / \tp(x + iy)\ 2 dy < +00 . 

x>a Jr. 

Applying the Laplace transform to both sides of the renewal equation 

f(t) = S(t)+[ K(s)f(t- s)ds, t>0 (18) 
Jo 

satisfied by /, we see that 

f(z)(l-K(z))=S(z), 5R(z)>0. (19) 

Since K £ L°°(R) with support in [0,2], the Laplace transform of K is entire 
holomorphic; it is given by 

k{z)= l - er2z{l 2 +2z \ z ^o, k(o) = i. 

2z z 

Since -^'(0) = — |, the only zero of 1 — K near z — is 0. Besides 

/•OO 

\K(z)\= e~ m[z) K{t)dt < 1 whenever $t(z) > , 
Jo 

since .ftT is a probability distribution on R + . Furthermore 

p 00 />oo 

Sft(A"(iy)) = / cos(ty)K(t)dt < / = 1 

Jo Jo 

so that z — is the only zero of 1 — K on the imaginary axis. Finally, for each 

a G R, one has 

IWI<^(l + (l + 2M)e- 23? «). (20) 

In particular | JiT(^r) | — > as 3(z) — > +00 uniformly in Si(z) > a. Therefore, there 
exists a > such that 

1 - k{z) = and ft(z) > -a <S> z = 0. 
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Hence z i— > - — ~—r is holomorphic on the half-plane II_ Q , and one recasts (ITT)1) as 

1 — K\z) 

zf(z) = zSi l) , Vt(z)>-a. 
1 — K (z) 

Indeed, since S £ L°°(R) with support in [0,2], its Laplace transform is an entire 
holomorphic function. 

Moreover, assuming that u m £ C*((0, 1)), we deduce from the formula 



S(t) = J u m (l +f- ts)K(s)ds 

that S £ C 1 (R+), with 

/OO 
(u m )'(l + t 2 - ts){2t - s)K{s)ds 

so that, for all i > 0, one has 

p OO 

\s'(t)\ < 2t\\(u m y\\ L ~ / K( s )d S < q(u in y\\ L o 



Since S(0) = u rn (l) = 0, one has zS(z) = S'(z); and since S' is continuous with 
compact support on R, its Laplace transform S' £ H 2 (H a ) for each a £ R. 

In view of (|2T)|) . for each a' < a, the function z i-> (1 — if (z)) -1 is bounded on 
n_ Q ' \ _B(0,a'/2), so that the function 

zS(z) 

z i y belongs to H 2 (IL- a >) 

1 — K(z) 

for each o[_< a. 

Since e at (j) — <j)(z — a), this implies the existence of a function g such that 
e a t g £ L 2 (R + ) for each a' < a and 

q(z) — — — , whenever ^R(z) > —a' . 
y ' l-K(z) V ' 



Consider the function G defined by 

G(i) 

Then G is continuous on R with support in R + , satisfies 



/ g(s)ds, i>0. 
Jo 



< \G(t)\ < 11 6 " , i>0, 
V2q;' 

so that G belongs to "H 2 (n^) for each 77 > 0. 
Now, for each 2 G II r; 

/•oo p oo pec 

zG{z) = I ze~ zt I g{s)dsdt = / g(s) / ze~ zt dtds 

JO Jo Jo J s 
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In other words, G and / are bounded continuous functions on R supported in R+, 
whose Laplace transforms coincide on ITo. Therefore 

f(t) = G(t) = [ g(s)ds , for each t > , 



■Jo 

and since e a l g G L 2 (R+) for each a' < a, one conclude that 

/ + oo / p+oa ^ V 2 

\g(s)\ds < We^gWv (J e^'ds 

- ^hle-^t < Ca ,e-<*'t 



V2a' 
for each t > 0. 

Going back to the original problem (fT5j) or (fl6|) . we deduce from the method of 
charateristics (fT7f that 

a;, = , whenever t > diam(f2) , 
or equivalently, in the frequency dependent case, 

diam(f2) 



I{t,x,u,v) = B v ({^f{t)) 1 ' 4 ) 



whenever t > 



since r x<ul < diam(J7). □ 

Before concluding, it is worth discussing the role of the monokinetic assumption 
in the exponential decay estimate obtained here. 

Intrinsically, the integral kernel K(s)ds in the proof of Theorem l5.1l is the image 
under oj r X)U of the probability measure proportional to (u> ■ n x ) + duj on S 2 . 
Since < t XiUJ < 2 = diam(fi), this kernel K is supported in [0,2]. Therefore 
its Laplace transform is holomorphic entire, and since solving the renewal integral 
equation is equivalent to the possibility of inverting 1 — K, the only constraint is the 
potential presence of zeros of 1 — K; besides, since K has integral 1, the half-plane 
3?(z) > is zero-free, except for z = 0, and zeros of 1 — K cannot accumulate on the 
imaginary axis. Hence the solution of the renewal equation is obtained by inverting 
the Laplace transform, deforming the contour of integration to be 3?(z) ~ —a' for 
some a 1 > 0. 

In the case of the collisionless gas in Theorem 14. 1[ the integral kernel K(s)ds is 
the image under v i— > t XiV of the probability measure proportional to (v-n x ) + M(v)dv 
on R 3 . Because of the low-speed gas molecules, its behavior near s = is the same 
as that of the image measure of, say, on the unit ball B(0, 1) under v t x , u 

— the Maxwellian weight M(v) obviously plays no role near v = 0. Thus, the decay 
of K for s — » +oo only depends on the Jacobian weight Iw^" 1 that appears when 
integrating radial functions with the Lebesgue measure in R N . This results in only 
algebraic decay for K in the limit s — > +oo, so that the Laplace transform of K 
does not have a holomorphic extension in any half-plane of the form 3?(z) > —a 
with a > 0. Therefore, the same deformation of contour as in the proof of Theorem 
15.11 is impossible in that case. In order to be able to push the integration contour 
in the Laplace inversion formula to the left of the imaginary axis, K should decay 
exponentially fast as s —> +oo. This would be the case if the Jacobian weight \v\ Ar_1 
appearing in the integration of radial functions with the Lebesgue measure in space 
dimension N could be replaced with e~ a /'"! for some a > 0. Thus, the difference 
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between the speed of approach to equilibrium observed in a collisionless gas with 
diffuse reflection of gas molecules on vessel walls kept at constant temperature, and 
the spectral gap property for the linear Boltzmann equation on the other hand, can 
be measured in terms of the vanishing order of \v \ N ~ X and e~ a /\ v][ as \v\ — > 0. 

6. Conclusion 

We have discussed in this paper the speed of approach to equilibrium for a colli- 
sionless gas enclosed in a vessel whose wall is kept at a uniform and constant tem- 
perature, assuming diffuse reflection on the vessel wall. We propose lower bounds 
for this convergence rate, depending on the possibility of initial distribution func- 
tions that are concentrated at low particle speeds. Not only do these lower bounds 
exclude exponential decay to equilibrium, they also constrain the power in any po- 
tential algebraic decay rate. On the other hand, assuming spherical symmetry of 
the initial velocity distribution function and of the vessel, we establish an algebraic 
decay rate by reducing the problem to a renewal equation for the flux of outgoing 
gas particles. All these results are consistent with the asymptotic behavior obtained 
in |10] with numerical simulations. To clarify the role of low-speed particles, we 
have studied the case of radiative transfer in the vacuum, in a container with infin- 
itely opaque boundary, so that the total incoming radiation flux on the boundary 
is reemitted as a Planck distribution. We establish exponential convergence of the 
radiative intensity to a Planck distribution in the long time limit, assuming the 
same spherical symmetries as in the case of a collisionless gas. 

There remain several open problems in this direction, such as a) establishing de- 
cay rates without any symmetry assumptions, and b) obtaining the optimal decay 
rates, even in the spherically symmetric case. Both problems might involve math- 
ematical techniques quite different from the ones used here. We hope to return to 
these questions in subsequent publications. 
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